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Abstract 

We present a simple model of a multi-quark droplet evolution based on 
the hydro dynamical description. This model includes collective expansion 
of the droplet, effects of the vacuum pressure and surface tension. The 
hadron emission from the droplet is described following Weisskopf's statis- 
tical model. We have considered evolution of baryon-free droplets which 
have different initial temperatures and expansion rates. As a typical trend 
we observe an oscillating behavior of the droplet radius superimposed with 
a gradual shrinkage due to the hadron emission. The characteristic life 
time of droplets with radii 1.5-2 fm are about 9-16 fm/c. 

1 Introduction 

In relativistic heavy ion collisions, a hot and dense fireball is produced in the 
overlap zone. It is believed that at early times this fireball is made of Quark-Gluon 
Plasma (QGP). Experiments aimed to study properties of QGP are going on at 
Brookhaven National Laboratory's Relativistic Heavy Ion Collider (RHIC) and at 
CERN's new Large Hadron Collider (LHC). Large amounts of data are already 
accumulated that helps to develop realistic theoretical models of the collision 
process and deconfinement-hadronization transition. No single approach has so 
far been able to describe all the aspects of this complicated process, however, 
several successful models have been proposed. Among the most popular is the 
hydrodynamical model, where the QGP is assumed to evolve as an almost perfect 
fluid. Such approach was first proposed by Landau almost 60 years ago p[], and 
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since that time its different versions are used to describe nuclear collisions at 
intermediate and high energies (see the recent paper [2] and references therein). 

One may expect, that in the course of fast expansion, this fluid will split into 
many droplets which later on evolve by evaporating hadrons from the surface [31 
Hj. The idea is that the strong collective expansion of matter, as seen in relativistic 
heavy-ion collisions, may lead to overcooling and subsequent fragmentation of 
the plasma phase into droplets surrounded by a low-density hadron gas. The 
dynamical implications of such a scenario were studied in refs. [S], El El El E] • In 
this paper, we will present a simple model of relativistic dynamics of a hot quark- 
antiquark droplet which may be created in this first step of the hadronization of a 
rapidly-expanding quark-gluon plasma. Our model is motivated by earlier works 
[TUl [TT] , where a simplified hydrodynamic model with hadron emission was used 
to describe the dynamics of a quark(-gluon) droplet. In our present approach we 
additionally consider effects of the surface tension and use the entropy balance 
equation for determination of the temperature. Also we present a detailed study 
of the droplet trajectories in collective coordinates. 

This paper is organized as follows. In Section 2, starting from hydrodynamic 
equations, we derive a simple model to describe evolution of an individual quark 
droplet. In Section 3, semi-analytical solutions of this model are obtained for 
the idealized case of no hadron emission. In Section 4, the hadron emission from 
the droplet is included following a statistical approach. In Section 5, we present 
results of numerical simulations of the droplet evolution taking into account the 
hadron emission. Conclusion and outlook are given in Section 6. 

2 Modeling evolution of a quark droplet 

2.1 Averaging hydrodynamic equations 

Let us consider a droplet consisting of deconfined matter at about critical 
temperature. As model calculation, show, under such conditions the gluons ac- 
quire a large mass and are strongly suppressed by the Boltzmann factor fTH [T5] . 
We take, therefore, only quarks and anti-quarks, in consideration. The quark- 
antiquark plasma inside the droplet is described by the hydrodynamic equations 
as a perfect fluid which is characterized by its rest frame energy density e, isotropic 
pressure P and spherically-symmetric velocity field v{r) with r being the distance 
to the center of the droplet. For the equation of state, we use, the MIT bag model 
i.e. an ideal gas of massless quarks and antiquarks confined in the spherical cavity 
of radius R(t). The energy momentum tensor of the fluid is represented as 

T'"' = Tj_a + TZnum = + " Pg'" + Bg^^, (1) 

where g'^'^ is the metric tensor g^" = diag[l, — 1, — 1, — 1] and = 7(1, {/(r)) 
is the collective 4-velocity. The last term in this expression corresponds to the 
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energy-momentum tensor of the "false" vacuum inside the droplet, characterized 
by the energy density B (bag constant). Let us first consider an idealized case 
where the droplet does not emit particles from the surface. Then the dynamics 
of the fluid follows the energy-momentum conservation equation for z/ = 0: 

a^T'^o = aoT°° + = 0, i = l,2, 3. (2) 

In the considered case the entropy is also conserved i.e. expressed by the equation: 

d^isu'') = 9o(sM°) + d^isu') = 0, (3) 

where s is the entropy density. Furthermore, we can add the continuity equa- 
tion expressing the conservation of quantum numbers such as baryon charge and 
strangeness: 

d^{nX) = ao(n,M°) + diiriju') = 0, j = B,S,... (4) 

where nj is the corresponding number density. From these equations, we will 
now derive the global dynamics of the droplet. Following ref. pUl [TT] this can 
be done by averaging hydrodynamic equations over the spatial coordinates. All 
our calculations are done in the center of mass coordinate system and under 
the assumption that the temperature and all thermodynamical functions are ho- 
mogenous within the droplet. For the velocity field of the fluid, we choose a 
Hubble-like profile of the form 

v = ^-R, (5) 

with R being the velocity of the surface of the droplet. A more general form of 
the flow profile as e.g. v = (|^) R (for positive a), could be considered too (see 
e.g. ref. [H]), but for the sake of simplicity, we will adopt the linear form. With 
this choice we average Eqs. ([2|), (p| ) and ^ over the volume V of the droplet. 
This yields for the l.h.s. of Eq. (|2j) 



{V[{e + P){^')-P + B]}. 



1 d 

vJt 



The second term in the r.h.s. of Eqs. ([6]) vanishes because the volume integration 
can be replaced by an integral over the remote surface, where all thermodynamic 
quantities vanish. Analogously, for Eqs. ^ and (|4]), we have 



1 /• „ _ I f d , ^ „^ I 

V J, 
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^^^d„(su'')dV^-J^-(s,)dV+-J^V-(sr^)dV 



(7) 



3 



and 



Eqs. ([g]), ([7]) and ([s]) contain spatial averages of 7 and 7^ (denoted as (7) and 
(7^)). Using the parametrization in Eq. (|5|, one can get the exphcit expressions: 

(7) = = ~ (arcsin(i?) - i^v^T^) , (9) 

(7^) = ^-j^ = 1^ (arctanh(i?) - R) . (10) 

At i? — 0, these expressions are reduced to 

(7)-i + ^i?^ + ^i?^ + ---, {i')^i + Ir' + Ir' + --- (11) 

Since the droplets are surrounded by a dilute hadron gas or even by the 
physical vacuum, there should be a surface energy associated with the change 
of the vacuum condensates inside the droplet. This means, we should add to 
Eq. ([g]) a surface term, which is commonly disregarded in the hydrodynamical 
calculations. In our present calculations we use a standard parametrization of 
the surface energy with form AirR'^a, where a is the surface tension coefficient. 
For the moving surface this expression has to be modified, because it should also 
contain the kinetic energy of the moving surface. The correct expression can be 
obtained by using an analogy with a relativistic particle of mass AirR^a moving 
with velocity i?, i.e. 

■ A-nR^a 

Esurfa.ce{R, R)= . (12) 
1 



Finally we obtain the following expressions for the total energy, entropy and any 
conserved quantum number in the droplet: 

E{R, R) = ^R' [ie + P) (7^) - P + B] + (13) 

SiR,R) = ^Rh{^), (14) 
Ayr 

N^{R,R) = —R\{^), J=B,S,... (15) 

Additionally, we need explicit expressions for internal energy density, pressure 
and entropy density for the plasma inside the droplet. For massless quarks and 
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antiquarks, these thermodynamic quantities as functions of temperature T and 
chemical potential /i are expressed as 



— I 

!0 



e(T, /i) = ^z.,T^ + -fT'^^' + ^/.^ P(T, /i) = -e(T, /.), (16) 



.(T,/.) = ^z.,T3 + ^T/.^ (17) 

nB(T,/.) = ^TV+^g, (18) 

where z/g is the degeneracy factor of quarks (antiquarks). Since we include: up, 
down and strange quarks, Ug = 18 in our calculations. For strange quarks Eqs. 



(16), (17) and (18) should be generalized to include finite rest mass ~ 150 
MeV. But for the sake of simplicity, in our calculations, we will disregard the 
strange quark mass. Please note, that in the thermodynamic expressions below, 
the variables T and fj, are often suppressed. 

2.2 Choice of model parameters 

To perform further calculations, we should specify the model parameters and 
initial conditions for the droplet evolution. As already stated above, for calcula- 
tion of the droplet bulk energy we use the MIT bag model. One should bear in 
mind, however, that in an equilibrated system this model predicts a first-order 
phase transition between the quark-gluon plasma and a hadron gas at all T and 
fi, see e.g. ref. |J^. On the other hand, the QCD based lattice calculations show 
a crossover type of the deconfinement phase transition at small fi. It occurs in 
the temperature interval (160-190) MeV, depending on the numerical scheme and 
quantity considered, for details see refs.p^ [T7j. In a non-equilibrium situation, 
when the transformation from quark-gluon to hadronic degrees of freedom is not 
fast enough, a formation of quark droplets is quite likely to happen. As shown 
in refs. [31 E], the QGP phase produced in relativistic heavy-ion collisions may 
split into droplets due to a strong collective expansion. Below, we assume that 
deconfinemed matter inside these droplets may exist in a metastable state at 
temperatures as low as 150 MeV and /x ~ 0. To describe this state we apply the 
bag model with bag constant B = 200 which is slightly lower than the stan- 
dard value around 235 With this choice, the balance between the plasma 



pressure (16) and the vacuum pressure (Pvac = —B) is achieved at a temperatures 
To ^ 150 MeV. 

For small droplets, the finite-size corrections to the droplet energy should be 
taken into account in addition to the bulk term proportional to R^. As shown in 
ref. [IS], within the MIT bag model with massless quarks, the standard surface 
term proportional to vanishes, and the curvature term proportional to R has 
a negative sign. The implications of this term for finite-size quark-gluon droplets 
in the mixed phase of the deconefinement phase transition, were studied in ref. 
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[T^ . The author's conclusion is that this correction is rather small (around 30 
%), even for droplet radii about 1 fm, and negligible for radii R^2 fm. 

The MIT bag model assumption of an idealized sharp boundary of the quark 
bag might not be realistic enough to correctly describe the surface properties. 
Other approaches, which are based on effective field-theorectical models, usu- 
ally predict a non-vashing surface energy. This energy is associated with the 
gradients of the mean fields (condensates) in the transition region between in- 
side and outside of the droplet. For instance, the calculations within the linear 
sigma model with quarks, ref. [20], predict the surface tension coefficient in the 



range (5 - 15) Other estimates using lattice QCD results |2I1 [22] give even 
larger values, of order of 100 Taking into account these uncertainties, in 

our exploratory study, we choose the surface energy in the form of Eq. (12) with 
temperature-independent coefficient o" = 50 

For the initial droplet radius we take the values of about 2 fm, which is 
motivated by estimates of ref. There the droplet size was related to the 

expansion rate of the quark-gluon plasma just before the break-up. It is clear, that 
this primordial plasma expansion should lead to the residual droplet expansion, 
hence in our investigation we also consider droplets with initial radial velocities, 
i.e. Ro ^ 0.2c 



3 Droplet dynamics without hadron emission 

3.1 Energy functional 

When hadron emission from the droplet is disregarded, its energy, entropy 
and net baryon number are conserved i.e. 



E{R, R) = Eq, 
S{R, R) = So, 
Nb{R,R) = Nb,o. 



(19) 
(20) 
(21) 



In the simplest case of zero chemical potential, we are left with expressions of the 



energy and entropy densities as functions of temperature only. Then Eq. (20) 



can be used to eliminate the temperature from Eq. ( 19 ) leading to an expression 



containing R and i? in a very nonlinear way. The resulting energy functional is 



^ ^' ^ = P o \4/3 + Q^^^ + / 

R 3(7)4/3 3 Vl-i?2 



(22) 



where (7) and (7^) are given by Eqs. (9) and (10), and k, 
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serting this expression in Eq. (19) leads to a nonlinear differential equation. 
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Figure 1: The equilibrium radius of the droplet as a function of temperature calcula- 
tions for a = 30,50 and 70 and B = 200 In our calculations, we primarily 
choose (indicated by the dot), cj = 50 S = 200 ^ and i?o = 2 fm which 
corresponds to a temperature of 153.4 MeV. 



whose direct solution is difficult to obtain. To solve it, we need further simplifi- 
cations, which are introduced below. First, let us analyze the shape of the energy 



functional (22). The equilibrium radius of the droplet, i?o, corresponds to the 



minimum energy state at i?, i.e. 



dE 



= —-^+ATTR'oB + 8TrRoa = 0. (23) 

R=Ro,R=0 ^0 



This condition is equivalent to the pressure balance at the surface 

^o(To) = B + ^, (24) 
where Pq is the pressure of the quark gas inside the droplet, B is the vacuum 



pressure, and ^ is the Laplace pressure. The numerical solution of Eq. (24) is 
presented in Fig. [T]for three different values of a. For the default values of Sand 
a, an equilibrium radius of 2 fm requires the initial temperature of 153.4 MeV. 
Then the minimal energy of the droplet is around 34.4 GeV. Obviously, to have 
higher initial temperature with the same radius, one should increase a or B. 

For further calculations, it is useful to make a graphical representation of the 
energy functional in phase space defined by variables R and R. The corresponding 
3-dimensional plot is given in Fig. [21 As one can see, the energy surface is 
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Figure 2: The droplet energy as a function of its radius and surface speed. Calculations 
are done for initial temperature: Tq = 154.3 MeV and model parameters a = 50 
and 5 = 200 The minimum corresponds to about 34.4 GeV at radius 2 fm. The 

colored contours corresponds to the iso-energetic contours at four different energies: 
34.6, 39.6, 44.6 and 54.6 GeV shown in the next figure. 



symmetric in the velocity coordinate R, while it is rather asymmetric in the R 
coordinate. The horizontal cuts of this surface at Eq = 34.6, 39.6, 44.6 and 54.6 
GeV are shown in Fig. |3j For lower energies, they look rather symmetric with 
respect to the point {Rq, 0) both in radius and velocity coordinates, while at 
higher energies, one can notice deviations at small radii. 

Rather surprisingly it turns out that these contours fit quite well with the 
shape of superellipses (see ref. [23]), defined by the equation 



R — Rc 


n 


R 


a 


+ 


~b 



(25) 



where Rc-, a, h and n are fitting parameters. These fits are shown together with 
the data points in Fig. [3} and the corresponding parameters are listed in Table 

m 

With parametrization of Eq. (25), we can easily find a numerical solution for 
the time-dependent radius of the droplet using standard methods like the Runge- 
Kutta scheme. In Fig. |4| the corresponding solutions are presented for the same 
total energies as in Fig. |3} The radial oscillations of droplets are explained as 
follows. Since we start with the same radius i?o = 2 fm and temperature Tq = 
153.4 MeV, the different total energies correspond to different initial velocities; 
around 0.14 for the droplet with smallest energy and 0.91 for the droplet with 
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Figure 3: Iso-energetic contours in {R, R) plane shown for four different energies (34.6, 



39.6, 44.6 and 54.6 GeV) together with fits to super ellipses (Eq. 25). The parameters 
for the fits can be found in Table [TJ 



the biggest energy (see Eq. 22). A droplet with an outward velocity will expand 



isentropically until the temperature of the thermal gas drops to a level where 
its pressure becomes smaller than the sum of the bag pressure and the Laplace 
pressure (see Eq. (24)). At this point, the droplet begins to contract until the 



plasma inside the droplet is heated up sufficiently to increase the thermal pressure. 
Generally, these solutions describe anharmonic oscillations of the droplet around 
the minimum energy state. The anharmonicity effects are getting more and more 
significant with increasing energy. The deviations from the harmonic oscillations 
are clearly seen in the shape of the curves in Fig. |4j The periods of the ocillilations 
are listed in Table [H 



Table f : Fitting parameters: a,b and n for the super ellipse equation together with 
the coordinate Rc at the center and the period of oscillation. 



E, [GeV] 


a, [fm] 


6, [c] 


n 


i?c, [fm] 




[c] 


r, [fm/c] 


34.6 


0.f4 


0.f4 


2.00 


2.00 


0.14 


6.19 


39.6 


0.65 


0.61 


2.10 


2.04 


0.61 


6.39 


44.6 


0.90 


0.77 


2.22 


2.08 


0.77 


6.68 


54.6 


f.21 


0.91 


2.65 


2.17 


0.91 


6.82 
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Figure 4: Numerical solutions for tlie radius of tlie droplet as a function of time 
obtained after fitting the iso-energetic contours to super ellipses. The coloring matches 
the previous plot. The curves labeled with numbers 1,2,3 and 4 correspond to energies 
32.5, 37.5, 42.5 and 52.5 GeV respectively. 



3.2 Small amplitude oscillations 



It is instructive to compare numerical results presented above with an an- 
alytical solution, which can be obtained in the case of small-amplitude oscilla- 



tions. Expanding the energy functional of Eq. (22) around the equilibrium point 
R = Rq, -R = 0, we get 



E{R, R) = — [3-3x + x^] 
Ro 

4 



1 + - {RoxY + ■ 
5 



+ -ttRIB [1 + 3x + 3x^ + ■■■] + ATiRla [l + 2a; + x^] 



(26) 



l + -(i?o±)' + 



where x = — 1^ Retaining only quadratic terms in x and x, we get 

E{R, R) = E{Ro, 0) + Cix^ + C2Rlx'^, (27) 

where Ci = + AitRqB + AnRla and C2 = + 2TTRla. In this approxima- 
tion the energy functional takes the form of a harmonic oscillator with angular 
frequency 



C2 Ro 




1 
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surface 



plasma 



K 



Rla, 



(28) 



where -^surface = 4:71 R^a and -Epiasma = Now we can estimate the oscillation 
period for the considered droplet (-^surface = 2.5 GeV, -Epiasma = 25.1 GeV ) as 

r = ^ = 6.11^. (29) 
cu c 



This fits well with our numerical solution for the lowest energy (Fig. [3]), which 
shows the period of 6.19 fm/c. At higher initial energies anharmonic effects 
become apparent and the oscillation period increases. 



4 Statistical description of hadron emission 

In this section, we consider the hadron emission from the droplet's surface. 
The idea is that the quarks close to the surface can combine into colorless hadrons, 
which may leave the droplet carrying away the energy Uh- We will describe this 
process following Weisskopf's statistical model in [24J. A similar approach has 
been used earlier in ref. [TT]. The ability to emit hadrons is determined by the 
thermal excitation energy E*[T, R, R) = E[T, R, R) — E{0, R, R). The droplet is 



assumed to be spherical at all times and in thermal equilibrium before and after 
each emission. According to Weisskopf's model, the double-differential emission 
rate is given by the ratio of the statistical weights before and after the emission. 




Figure 5: Hadron emission from the surface of a quark droplet of radius R. 
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dp dt 2 27r2 



1 2 

r 

2 47r2 



pMexp[5(E^ iV, i?) - iVo, i?o)] 

(30) 



where p is momentum of emitted hadron, Vh is its spin-isospin degeneracy factor, 
A is often estimated as the cross section of inverse reaction, but for our problem 
it is more appropriate to identify it with the total surface area of the droplet, A = 
-^geom = 47ri?2, S{E'^,R) is the entropy of the droplet at the thermal excitation 
energy E* and volume V = ^nR^. Note that we have introduced a factor of 1/2, 
since the particle emission is only possible for polar angles ^ < f • Furthermore, 
in the case of a moving surface, the emission rate will be suppressed by the 7 
factor ^llj . The energy conservation in the emission process can be expressed as: 



AEh = E-Eo = -Uh, Uh = ^Jml+p\ (31) 

where m/j is the hadron mass. The change in the entropy due to emission of a 
hadron h can be found from the 2nd law of thermodynamics: 

A^^ = ^ [AEl + PAVh - fih^N,,) , (32) 

where AE^ = —ooh, P is the thermal pressure and AV^ is the change of the volume 
due to emission of one hadron, which will be neglected below, ANh = — 1 since one 
particle is emitted and is the chemical potential of the droplet corresponding 
to the quantum numbers of the emitted hadron. Now, the double-differential 
emission rate can be written as 



+ 



T T 



(33) 



The first two moments of this distribution give the total particle number and 
energy emission rates 



dt vr 



i?^ / p^e T dp, (34) 







— y ^ — R I p ym^+ p^e t dp. (35) 



dt ^ vr 

h 







In our calculations, we include emission of mesons {tt, K, p,u), baryons (A^, A) 
and hyperons {Y = A, S). The Fermi and Bose statistics effects are disregarded 
in these calculations. Obviously, the energy loss due to emitted particles leads 
to a change of the entropy and total energy of the droplet. Therefore, instead of 
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Eqs. (19) and (20) we should now solve the following equations 



dE 
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~dt ~ 


dt 


dS 


d 


~dt ~ 


dt 



V{{e + P){j')-P + B} 



Vl -i?2 



dE^ 



loss 



dt 



[Vs{j)] 



1 dE, 



loss 



T dt 



(36) 
(37) 



These equations together with Eq. (35) should be solved for functions R{t) and 
T{t). 



5 Numerical simulations of droplet dynamics 



To solve Eq. (35), (36) and (37), we have created a numerical scheme calcu- 
lating at each time step (n ■ At) the change of the variables R, R and T based 
on their previous values. We are using the Runge-Kutta method to determine 
Rn+i and Rn+i from Rn, Rn and Rn- The acceleration Rn can be calculated by 
knowing S'„, Rn and Rn- Now T„+i can be calculated from S'„, Rn+i and Rn+i- 
Finally, from calculating the energy loss from Tn+i and Rn+i, we get Sn+i- Now 
the loop can be redone at the next time step: [n + 1) ■ At. We have performed 
numerical simulations for four initial conditions, all with -Ro= 2 fm : 1) -Ro=0, Tq 
= 153.4 MeV (close-to-equilibrium case), 2) Ro=0.2, Tq = 153.4 MeV, 3) i?o=0. 
To = 140 MeV and 4) Ro=0, Tq = 175 MeV. 



5.1 Evolution of droplets with hadron emission 

The behavior of the droplets with hadron emission differ significantly from 
that in the emission-free case. In Fig. [6} the surface velocity is plotted as a 
function of radius of the droplet, analogous to Fig. |4]for the emission-free case. 
The dynamics shows a rather peculiar oscillating behavior superposed with the 
shrinkage of the droplet. For the close-to-equilibrium case without initial speed, 
the droplet does not acquire any outward speed at all. In other words, the general 
trend is a shrinkage of the droplet due to the hadron emission. Nevertheless, in the 
case of high initial temperature or non-zero initial speed, the droplet is initially 
expanding until -R = at radii of about 2.32 fm and 2.16 fm, respectively, and 
then shinks again. 

Fig. [7[ shows the radius of the droplet as a function of time. The solid line, 
corresponding to the close-to-equilibrium initial state shows, a gradual decrease of 
the radius due to the hadron emission. When the initial state is out of equilibrium 
(3 other curves), the droplet 's radius exhibits damped oscillations, reminiscent of 
the droplet oscillations in the emission- free case (see Fig. |4]). In all cases, after 
averaging over oscillations, the radius drops almost linearly with a rate of about 
0.075 c. We stop our calculations at radii of about 0.8 fm, close to the charge 
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Figure 6: Trajectories of tlie quarlt droplets in tlie R — R plane are shown for four 
different initial condition (indicated in the figure). A cut off at i? = 0.8 fm is applied. 
The dots and the arrows on the plot indicates the initial points and the direction of 
the droplet evolution. 

radius of the proton, because at smaller values our macroscopic considerations 
cannot be justified. As one can see, from Fig. [7} a droplet with initial radius of 2 
fm shrinks to this size in about 16 fm/c. This is of the order of the hadronization 
times extracted from the HBT measurements, which are about 10 fm/c at LHC 
and around 7 fm/c at RHIC |2Sj- However, one should bear in mind that the life 
time of the droplet is very sensitive to its initial radius. For instance, a droplet 
with radius of 1.5 fm will hadronize within the time interval of about 9 fm/c. 

In Fig. [8| the temperature of the quark gas inside the droplet is plotted as 
a function of time for the same four initial conditions as discussed above. One 
can see, that for the high-temperature case and the non-zero initial speed case 
the temperature of the quark gas decreases and reaches around 140 MeV at the 
turning point for both cases. At that point, the thermal pressure is too low to 
balance the vacuum pressure and the Laplace pressure, and the droplet therefore 
starts to contract until the temperature reaches a new maximum around 175 MeV 
for both cases. 

We can make two conclusions from Fig. [8] First, the temperature of the 
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Figure 7: The droplet radius as a function of time for four different initial conditions 
(indicated in the figure). A cutoff of R=0.8 fm is applied. 



droplet slowly rises during its life time, and second, oscillations of the tempera- 
ture get smaller with time. This shows that due to the damping, the temperature 
needs to increase slightly to balance the surface tension of the smaller droplet i.e. 
the droplets are evolving towards an equilibrium state, which however, is con- 
stantly changing to higher temperatures according to the curve depicted in Fig. 
[l] This trend is different from the one found in ref. [llj, where the surface tension 
was disregarded, so that the droplet temperature was approximately constant. 



5.2 Characteristics of emitted hadrons 

Now, let us discuss characteristics of the emission process. We start with 
the energy budget of the droplet. Fig. [9] shows, the total energy of the droplet, 
the thermal energy of the quark gas, the surface energy and the vacuum energy 
as functions of time for two different initial conditions i.e Rq = (solid lines) 
and i?o = 0.2 (dashed lines) at the same temperature To =153.4 MeV. The two 
cases differs only by around 1 GeV, since the difference in initial speed is rather 
small. At all times, the dashed line lies above the solid line for the total energy 
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Figure 8: The droplet temperature as a function of time for four different initial 
conditions (indicated in the figure). A cutoff at R=0.8 fm is applied. 



as we would expect. In the case of non-zero initial speed, we see a more violent 
dynamics in accordance with the earlier plots. In both cases, we have an initial 
energy of about 40 GeV which drops to around 4 GeV during the lifetime of the 
droplet. Most of the energy is accumulated in the thermal gas; about 30 GeV 
out of 40 GeV at initial stage. One can also see, that the oscillations of the total 
energy are coupled to the variations of the thermal energy. The oscillations of the 
vacuum energy are entirely due to the oscillations of the droplet volume, while 
the oscillations of the surface energy are caused by variations of the surface area 
and the speed. 

Figures 10 [11 and 12 show the hadron emission rates, energy loss rates and the 
number of hadrons emitted as functions of time for two different initial conditions 



as in Fig. M While in Fig. 12, only the zero- velocity case is shown. From Figs. 



10 and ^2] one can see that, the pious and kaons, are by far the most abundant 



emitted particles, since the emission of heavier species is suppressed by their mass. 
Nevertheless, the heavy particle spieces contribute significantly to the energy loss. 



It follows from Figs. 10 and n^Tl that at initial stage the typical emission rates 



are about 3 pious, 1 kaon and 1 heavier hadron per fm/c. These particles carry 
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Figure 9: Time evolution of different contributions to the total energy of the droplet, 
the thermal energy of the quark gas, the surface energy and the vacuum energy due to 
the bag constant. The solid and dashed lines represent initial conditions with Rq = 
and Rq = 0.2 respectively. 



away accordingly 1 GeV, 0.8 GeV and 0.8 GeV of energy per fm/c. From Fig. 
12 one can see that about 30 pions, 10 kaons and 7 heavier particles are emitted 
from the droplet with initial thermal energy of 30 GeV during its life time. 

From this analysis, we can make two conclusions: First, the inclusion of heav- 
ier particles in the emission process leads to a significant enhancement of the 
energy emission rate, and second, even moderate initial collective velocity leads 
to strong droplet oscillations and pulsed hadron emission. 



6 Conclusion and outlook 

We have presented a simple hydro-kinetic model of a quark droplet evolu- 
tion, which takes into account bulk and surface contributions to the energy, col- 
lective expansion and hadron emission from the surface. If hadron emission is 
disregarded, the droplet behavior is described by anharmonic oscillations. We 
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Figure 10: Emission rates of different hadron species (indicated on the figure) and the 
total emission rate as function of time for two different initial conditions. The solid 
and dashed lines represent calculation with Rq = and Rq = 0.2 respectively. 

have found an approximate solution for these oscillations using a super-ellipse 
parametrization of the phase-space trajectories. In the limit of small-amplitude, 
the period of oscillations is also proportional to the droplet radius. When hadron 
emission is included, the motion is changed to damped oscillations characterized 
by the decay time, which in first approximation is proportional to the droplet 
radius. For droplets with initial radii 1.5-2 fm these lines lie in the interval 9-16 
fm/c. Our calculations show that pions account for about 2/3 of the emitted 
hadrons, but they carry away only 50 % of the energy. 

In this paper we have considered only baryon-free droplet made of quarks 
and antiquarks with zero chemical potentials. This is a reasonable assumption 
for plasma produced at RHIC and LHC energies |2S]- At lower collision en- 
ergies effects of baryon asymmetry become more and more important. So, in 
the future, we are planning to include non-zero chemical potentials and investi- 
gate asymmetry between baryon and anti-baryon emission. Then, we can study 
such interesting phenomena as cold stranglet formation and flavor distillation 
[m |271 [2H] • The possibility of instantaneous hadronization of the dilute plasma 
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Figure 11: Energy emission rates for different hadron species (indicated on the figme) 
and the total energy emission rate as function of time for two different initial conditions. 
The solid and dashed lines represent calculation with Rq = and Rq = 0.2 respectively. 



[1] should be considered as well. 

In our calculations, we have made an assumption that the droplet keeps a 
spherical shape at all times. This assumption can of course be questioned because 
the recoil effect from the emission, the Rayleigh- Taylor instability due to viscosity 
and collective quadrupole modes could lead to an asymmetric shape or even a 
bulged or fingered surface of the droplet. The recoil effect due to hadron emission 
can most hkely be dismissed due to a stochastic nature of the emission process, 
i.e. the random directions of the emitted particles will lead to a cancellation of 
the recoil momenta due to individual particle emissions. Also, the surface tension 
should work against a non-spherical shape, trying to minimize the surface area 
of the droplet. 

For a more realistic description, one could introduce viscous terms in the 
energy-momentum tensor, and study their impact on the dynamics of the droplet. 
It is clear, however, that these terms will damp collective oscillations but not 
change dramatically the life time of the droplets, which is mainly determined by 
the hadron emission rate. Also, the finite-size corrections to the droplet energy 
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Figure 12: Total numbers of specific liadrons (indicated in the figure) emitted as a 
function of time. The initial conditions are To = 154.3 MeV, Rq = 2 im and Rq = 0. 

should be studied in more detail, in particular, higher-order terms in ex- 
pansion should be introduced. These terms could play an important role at late 
states of the droplet evolution. 
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